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Octave 


► http://www.gnu.org/software/octave/ 

► High-level language for numerical computations 

► Mostly compatible with MATLAB 
(http://www.mathworks.com/) 

Just about everyone thinks that the name 
Octave has something to do with music, but it is 
actually the name of one of the author's former 
professors who wrote a famous textbook on 
chemical reaction engineering, and who was 
also well known for his ability to do quick "back 
of the envelope" calculations. We hope that this 
software will make it possible for many people 
to do more ambitious computations just as 
easily. 
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Constant-Q Transform 


► [Brown, 1991] 

► One way to get a log-frequency spectrum 

► 0 = f/Af, frequency over bandwidth 

► DFT has constant bandwidth, so Q changes with 
frequency 

► Constant-0 transform varies bandwidth with 
frequency 

► Bandwidth is inversely proportional to number of 
samples 


Constant-Q Transform Definition 


N k cq — 1 

Xcq(ftcq) = £ w t Jn) X (n)e- l2 ' 0n,N ‘« 

fccc i n =0 

► 0 = l/(2 1 / b - 1) with b bins per octave (12 for 
semitone spacing) 

► f kcq = 2 kc ^/ b f m - m is the center frequency of the Ac cq th 
bin. 

► N kcq = Qf s /f kzq is the number of samples 

► w k is a window function (e.g. Hanning) 



Brown-Puckette Algorithm 


► [Brown and Puckette, 1992] 

► Can be more efficient to compute 

► FFT x (sparse) spectral kernel 

N 0 — l 

^cq(kcq) = jj- ^ X(k)K(k C q,k) 
0 k =0 
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Spectral Kernel 


N o — 1 

K(k cq ,k) = J2 w kcq {n)e' 2vQ/Nk ^e' 


n =0 
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Code 


Download code at 

http://hans.fugal.net/research/cq-octave/ 
cq Original 
ecq Brown-Puckette 
vecq Vectorized Brown-Puckette 
secq Brown-Puckette with sparse kernel 

svecq Vectorized Brown-Puckette with sparse 
kernel 

cqkern Spectral kernel 

sparse Convert spectral kernel to sparse matrix 


cq 


N kcq -i 

Xcq(*c) = E "tc q (n)x(n)e-' 2 '°" /w >'« 

fcc( i r)=0 

1 for kcq = 0:B-1 

2 f = f_min * r^kcq; 

3 Nkcq = round(Q*sr/f); 

4 Xcq(kcq+1) = sum (hamming(Nkcq) .* (x(l:Nkcq) \ 

5 .* exp( j2pi0n(l:Nkcq)/Nkcq))) / Nkcq; 

6 end 
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ecq 


1 X = fft(x, 

2 for kcq = 

3 Xcq(kcq) 

4 end 


N o-l 

^cq(^cq) = ^ X(k)K(k cat ,l<) 

0 k =0 

N0); 

1: B 

= sum(X .* K(kcq, :)) / N0; 
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vecq 


^cq(^cq) 


1 
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E mK(k cai ,k) 

k= 0 


1 X = fft(x, N0); 

2 Xcq = sum((repmat(X, B, 1) .* K), 2) / N0 


Parameters—First Experiment 


► f s = 44100 Hz 

► b = 12 bins per octave (semitone spacing) 

► f min = 16.35 Hz (CO) 

► (max = 22050 Hz (Nyquist frequency) 


Parameters—Second Experiment 


► f s = 44100 Hz 

► b = 12 bins per octave (semitone spacing) 

► f min = 130.81Hz (C3) 

► (max = 22050 Hz (Nyquist frequency) 


Specifications 


MacBook 

► 2 GHz Intel Core Duo 

► 2 GB RAM 

► OS X 10.5 (Leopard) 

► Octave 3.0.3 from MacPorts 
Athlon 

► Athlon 64 2800+ (1.8 GHz) 

► 1GB RAM 

► 32-bit Ubuntu 8.10 

► Octave 3.0.1 stock 
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Results, f min = CO 



cqkern sparse 


► cq is fastest 

► Vectorized are slower 

► Sparse is faster 
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Results, f min = C3 



► cq is slowest 

► Vectorized depends on hardware 

► Sparse depends on hardware 

















Discussion 


► Brown-Puckette algorithm not always more efficient 

► Same order: O(/V 0 log/V 0 ) 

► More space: O(N 0 log N 0 ) vs. O(N 0 ) 

► Vectorization worse or marginal, probably not worth 
it 

► Sparse marginal improvement, may be worth it 


Future Work 


► MATLAB 

► relative and absolute comparison 

► repmat 

► C 

► fast loops 

► more control over memory 
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